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ABSTRACT 

We studied superclusters of galaxies in a volume-limited sample extracted from the 
Sloan Digital Sky Survey Data Release 7 (SDSS/DR7) and from mock catalogues based 
on a semi-analytical model of galaxy evolution in the Millenium Simulation. A density 
field method was applied to a sample of galaxies brighter than M r = — 21 + 51og/iioo 
to identify superclusters, taking into account selection and boundary effects. In order 
to evaluate the influence of threshold density, we have chosen two thresholds: the first 
maximizes the number of objects (Dl), and the second constrains the maximum su- 
percluster size to ~120 h~ 1 Mpc (D2). We have performed a morphological analysis, 
using Minkowski Functionals, based on a parameter which increases monotonically 
from filaments to pancakes. An anti-correlation was found between supercluster rich- 
ness (and total luminosity or size) and the morphological parameter, indicating that 
filamentary structures tend to be richer, larger and more luminous than pancakes in 
both observed and mock catalogues. We have also used the mock samples to compare 
supercluster morphologies identified in position and velocity spaces, concluding that 
our morphological classification is not biased by the peculiar velocities. Monte Carlo 
simulations designed to investigate the reliability of our results with respect to random 
fluctuations show that these results are robust. Our analysis indicates that filaments 
and pancakes present different luminosity and size distributions. 
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1 INTRODUCTION 

It is well known that galaxies are not randomly distributed 
in the Universe, with high-density regions being observed 
as large-scale structures and low-density regions as voids. 
Under the current ACDM cosmological paradigm, the evo- 
lution of these structures started in the early Universe from 
primordial density fluctuations just after inflation, leading 
to the observed cosmic web. At very large scales, of tens 
of Mpc, clusters, groups, and even pairs or isolated galaxies 
are disposed in very large associations, sometimes of filamen- 
tary or planar structure. These associations are the largest 
non-virialized structures in the Universe: superclusters of 
galaxies. Their dynamical future is still uncertain, but in a 
dark energy dominated Universe most of them may evolve 
to island univ erses, single, isolated an d highly concentrated 
mass clumps (|Arava-Melo et al1l2009l ). 

The study of very large scale structures started with 
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Ide Vaucouleursl l| 19531 ). who identified a high-density region 
in the galaxy distri bution on the sky, nowadays known as the 
Local Supercluster. lAbelJ l|l958l ) also helped to unveil large- 
scales through his catalogue of clusters of galaxies identified 
in the Palomar Observatory photographic plates; he defined 
superclusters as clusters of clusters of galaxies. Larger sur- 
veys were carried out and, as a consequence, the distribu- 
tion of galaxies in the local Universe could be studied in 
detail. The Harvard Center for Astrophysics (CfA) survey 
meas ured the redshift of a sample of galaxies brighter than 
14.5 (|Huchra et al.| [l983). showing the filamentary distribu- 
tion of galaxies, with galaxy clusters at the connection of 
the filaments. These red shifts allowed to constrain the cos- 
molog ical model. Indeed, Efstathio u. Sutherland fc Maddoxl 
(1990) showed, from the analysis of large-scale galaxy clus- 
tering in the IRAS survey, that a cosmological constant was 
required to explain the galaxy distribution in the framework 
of the CDM model. 

The data derived from redshift surveys later allowed to 
study t he propertie s of individual s uperclusters, like Pi sces- 
Cetus (|Tullvl 1 19881 ) and Shapley (|Proust et al.1 l200rj) . as 
well as those of the whole population of superclusters, re- 
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vealing that they tend to have an elo ngated morphology 
and extensions up to ~100 h _1 Mpc (jBahcall fc Soneiral 
1 19841). Recent works reveal extensions up to 110-130h _1 Mpc 
(|Pandev et al. 2010). The study of superclusters highly ben- 
efited from the 2 degree Field Galaxy R edshift Survey (2dF- 
GRS) (|Colless. Dalton fc Maddaxll200ll) and th e Sloan Dig- 
ital Sky Survey (SPSS) (|Abazaiian et al.ll2009h . 

Indeed, more complete studies of large-scale struc- 
tures were possible using these large redshift surveys. Jaan 
Einasto's group used the 2dFGRS data to g enerate a 
catalogue of superclusters (|Einasto et al.l l2007al hereafter 
E07a) . Th ey also compared obse rved superclusters to simu- 
lated ones l|Einasto et al.ll2007bl . hereafter E07b) and stud- 
ied the spectral prop erties of galaxies within superclusters 
l|Einasto et al.ll2007d ). Their main results indicate that the 
overall properties of simulated and observed superclusters 
present good agreement with each other, but their lumi- 
nosity and multiplicity (number of galaxies) distributions 
seem to be different. Additionally, they found that galaxy 
morphology in superclusters depends on their richness, with 
rich superclusters presenting an early-type fraction slightly 
higher than poor superclusters. In another series of papers 
(Einasto et al. 2007d,e), these authors have studied the 
richest superclusters identified in the observations. Compar- 
ing the dumpiness of simulated and observed superclusters, 
they conclude that the dumpiness of galaxies in simula - 
tions is different from that observed (jEinasto et al.ll2007dh . 
and that the global and local environments are quite im- 
portant for galaxy m orphology and star formation activity 
l|Einasto et al.ll2007el ). 

Since superclusters are non-virialized stru c tures , 
they present a v ariety of morphologies (|Westl 1 19891 ; 



IPlionis et al. 19921 ). Several studies have used shape 
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Shandarin 1998) and 



statistics l|Sahni. S athvaprakash 

Minkowski Functionals ijMecke. Buchert fe Waenerll 19941 ) to 
determine topological and geometrical properties useful for 
morphological analysis. The SDSS and PSCz (Saunders et 
al. 2000) superclust ers seem to hav e a prevalence of fila- 
mentary structures (|Basilakos|[2003h . as well as a concor- 
dance i_wit4i_the_ACDM_jnodel_^ forma- 
tion (jBasilakos. Plionis fc Rowan-Robinsonl l2001) (hereafter 
B01). A comparison of observed and simulated superclus- 
ters showed that simulated superclusters are very similar 
to those observed, but the number-density of very luminous 
superclus ters seems to be hig her in observations than in sim- 
ulations l)Einasto et al ] |2006l . hereafter E06). 

Morphological studies suggest that galaxies are found 
in two distinct classes of structures at very large scales: 
filaments and pancakes. Using the shap efinder technique 
l|Sahni. Sathvaprakash fc Shandarinlll998l ). B01 applied this 
approach to distinguish between these two classes using 
the so called shape-spectrum. Further works used fila- 
ment features to constrain the galaxy clustering, since 

the bias parameter is also sensitive to filamentarity 

ijBharadwai fc Pandevll2004h. Using g alaxy luminosities and 
colours, Pandev fc Bharadwail (|2006l ) found a dependence 
between galaxy properties and filamentarity, proposing a 
scenario where elliptical galaxies are predominantly in dense 
regions, while spiral galaxies are distributed along filaments. 
A strong spatial alignment between clusters and host super- 
clusters in large filaments was found in N- body simulations 
l|Basilakos et al.ll2006l : iLee fc Edvardll2007f ). 



Here we present a study of supercluster morphologies 
through the study of the galaxy distribution in volume- 
limited samples extra cted from SDSS Data Release 7 
l|Abazaiian et al]|2009T ). We use a kernel-based density field 
method to identify the superclusters and Minkowski Func- 
tionals to quantify their shape. 

This paper is organized as follows. Section [2] presents 
the data used here as well as our method to deal with selec- 
tion effects. In section[3]we describe the kernel-based density 
field method used to identify superclusters as well as the cri- 
teria to classify enhanced regions as superclusters taking into 
account the selection and boundary effects. In section [4] the 
morphological classification is described and in section [S] we 
discuss the morphology of observed and simulated superclus- 
ters. Finally, in section[6]we summarize the main conclusions 
of this paper. In Appendix [X] we discuss the sensitivity of 
our supercluster identification and morphological analysis to 
the adopted kernel. 

When necessary, distances were calculated assuming the 
following cosmology: Q, m — 0.3, S1a = 0.7 and Hubble pa- 
rameter iio = 100 h^QQ km s _1 Mpc _1 . 



2 DATA 

The analysis presented in this paper is based on a volume- 
limited galaxy sample extracted from the Sloan Digital Sky 
Surv ey Data Release Seven (SDSS-DR7) (|Abazaiian et all 
2009). We have considered galaxies with measured radial ve- 
locities and with absolute magnitudes in the r band brighter 
than —21 + 5 log ft in the redshift range 0.04 < z < 0.155. 
Absolute magnitudes were calculated with k-corrections ob- 
tained with the code KCORRE CT vA.lA and with a spe - 
cific SDSS package provided by iBlanton" fc Roweij l|2007fl . 
Since superclusters may extend over several degrees on the 
sky, we have considered only galaxies within stripes 10 to 
37 to assure a large continuous area on the sky. The total 
number of galaxies selected is 120,013. 

In Section 5 we shall compare some of our results 
with numerical simulations. For this we have u sed simu- 
lated light-cones produced bv lCroton et al.l (|2006l ). based on 
a semi-analytic galaxy evolution model applied to the out- 
put of the Millenium Simulation (|Springel et all 120051 ). We 
have selected the four light-cones with parameters suited 
for SDSS: SDSS.SAcone_012_000, SDSS_SAcone_012_100, 
SDSS_SAcone_120_000, and SDSS_SAcone_201_000. Each 
covers an area of 60 x 30 deg 2 and the simulated galaxies 
were selected following the same criteria adopted in the se- 
lection of our volume-limited sample of SDSS galaxies. The 
number of simulated galaxies selected in the four light-cones 
is 99,850. 



3 THE DENSITY FIELD METHOD 

Superclusters are sometimes defined as l arge-scale over- 
densi ty regions in the galaxy distribution (|de Vaucouleursl 
1953). Adopting this definition, the density field method rep- 
resents a convenient way to identify these structures (e.g., 
E07a, B01). In this section we describe how we define the 
density field of a sample of galaxies. 
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3.1 The density field 

Firstly, using the equatorial coordinates (a, 8) and the red- 
shift z of each galaxy, we calculated its cartesian coordinates 



x = d c cos (<$) cos (a) 
y — d c cos(<5) sin(a) 

z = d c sin(<5) (1) 

where d c (z) is comoving distance of the galaxy. 

The luminosity density of the galaxy distribution, D(r), 
is calculated through the kernel approach. At a certain point 
r in space it is given by 



D(p) = 2^iC(|r-ri|,<7)L i W i (r,) 



(2) 



where K(r, a) is the kernel used to smooth the galaxy dis- 
tribution, Li is the luminosity of the i-th galaxy (at position 
ri) and Wi(ri) is a weight which takes into account selection 
effects (discussed in the next section). 

We adopt here Epanechnikov's kernel, which mini- 
mizes the asymp totic mean integrated squared error (e.g. 
ISilverm ann 1986) and is defined as 



K{r,a) 



f[l-(>» 2 ] 







r < a 
r > a 



(3) 



Here a is the smoothing parameter. We have adopted 
a=8 h _1 Mpc. The reason is that the number density of 
galaxies in our sample (see next section) is n = 2.1 x 
10 _3 (h _1 Mpc) -3 , corresponding to a mean distance be- 
tween galaxies of ~8 h Mpc. As shown later, this choice 
leads to a density field relatively insensitive to peculiar ve- 
locities. The density field is sampled in a 3D grid with cells 
of side l ce ii—4: h _1 Mpc. 

To identify structures in the density field, it is neces- 
sary to define a density threshold to separate high-density 
regions (e.g., superclusters) from low-density regions (e.g., 
voids). In this way, we have rejected all grid points below 
the threshold. Afterwards, a friends-of-friends algorithm was 
used to connect nearby high-density grid points, assigning 
them to single objects. The linking-length used is equal to 
the diagonal of the cell grid, i.e., lf a f — v3 l ce u —7 h Mpc. 
Only objects with more than 10 galaxies and volume larger 
than two grid cells, V min =2(l C eii ) 3 =128(h~ 1 Mpc) 3 , will be 
considered in the analysis |E07a]). 

Note that about 0.3% of the DR7 imaging footprint area 
are marked as holes. In these regions we have used bilinear 
interpolation to obtain the density field, considering only 
grid points at a given redshift more distant that 8 h _1 Mpc 
from the hole borders. Since the area occupied by the holes 
is small, and most of the holes are in low-density regions, it 
can be verified that this procedure has a negligible impact 
on our results. 

There is no natural value for the threshold density. In 
Figure [T] we present the number of structures (hereafter 
called superclusters) as a function of the threshold in units 
of mean density (£>), computed with the selection function 
discussed in 3.2. For low threshold density values, percola- 
tion links distinct structures and consequently the number 
of superclusters is low. At high threshold values, only high 
density objects are identified, also resulting in a low number 




D/D 



Figure 1. Number of superclusters as a function of the threshold 
density (in units of the mean density). 



of structures. In this work we have adopted two distinct val- 
ues for the threshold density. The first one, Dthresh — 3 x D 
(hereafter Dl), is the value which maximizes the number of 
superclusters. The second one, Dthresh = 6 x D (hereafter 
D2), was chosen such as the largest superclusters_present 
an extension of ~120 h _1 Mpc, as adopted by E07a. This 
length consists of the diagonal of the b ox which contains a ll 
galaxies of the supercluster, i.e., I — y Ax 2 + Ay 2 + Az 2 . 

Two important features of superclusters can be defined 
here: their richness and total luminosity. The richness can 
be written as 



R= Y^Wi 



(4) 



where N ga i represents the number of galaxies of the super- 
cluster and Wi is the selection effect correction of the i-th 
galaxy. The total supercluster luminosity- actually the ex- 
pected luminosity above the magnitude limit- is defined as 



Ltot 



"gal 



(5) 



where Li represents the luminosity of the i-th galaxy. 

It is worth mentioning that our results are not too sen- 
sitive to the choice of the smoothing kernel. We present in 
Appendix|X]a summary of results obtained with a truncated 
Gaussian kernel, showing that our estimates of supercluster 
parameters are indeed very robust. 

3.2 A model for the selection function 

The selection function aims at correcting for galaxies 
brighter than our magnitude limit that, for a reason or an- 
other, were not included in the sample. 

Indeed, our magnitude limited sample is affected by in- 
completeness due to fiber collisions in the spectroscopic sur- 
vey. Consequently, although the nominal magnitude limit of 
the SDSS Main Galaxy Sample MGS) is m r = 17.77, not 
all galaxies brighter than this limit were observed. There is 
a minimum distance between fiber allocations by the SPSS 
spectrographs of about 55 arcsec (|Strauss et al] 12002) and 
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Figure 2. Fraction of galaxies with observed spectra as a function 
of apparent r-band magnitude in the range 13.0 < m r < 17.77. 
The red line represents the fourth order polynomial fitted to the 
observed trend. 

some galaxies within the MGS photometric limits do not 
have spectroscopy because they are closer than 55 arcsec to 
a galaxy to which a fiber was allocated. This spectroscopic 
incompleteness depends on the apparent magnitude since 
fainter galaxies are more affected. This is shown in Figure [2l 
which presents the fraction of galaxies with spectroscopy as 
a function of the apparent magnitude m r . 

The spectroscopic incompleteness leads to a radial se- 
lection effect, since galaxies with higher apparent magni- 
tude tend to be at higher redshifts. Since in this case 
the shot-noise increases with distance, the coupling be- 
tween this radial effect and a constant smoothing param- 
eter introduces an additional distortion in the d ensity field 
|Gaztanaga fc YokovamaFl 993; Scli ak et al . 2009). This bias 
leads to an overestimation of the density with increasing red- 
shift (B01). 

To deal with these effects, we have adopted a sim- 
ple model for the selection function, with two components. 
The first one depends on the apparent magnitude (Si(m r )) 
and the second one on the redshift, or comoving distance 
(S2(d c )). The selection function is thus defined as 

S(m r ,d c ) = S 1 (m r )S2(d c ) (6) 

and is related to the weight W as 

w = sk.4)" 1 

The apparent magnitude component can be denned as 
Si(m r ) = n 3pec /ntot and we model the trend seen in Fig- 
ure [2] with a fourth order polynomial: 

Si(m r ) = 0.588605 - 1.941834m r + 0.419142m^ 

-0.029956m? + 0.000724™* (7) 

To estimate the radial component of the selection func- 
tion we, initially, calculated the mean of the density field grid 
points in ten regions with same volume, taking into account 
only the apparent magnitude incompleteness (i.e., assuming 
S2(d c ) — 1). As shown in Figure the mean density of each 
region increases with redshift, reflecting the bias mentioned 
above. However, a comparison with Figure 3 of B01 indi- 



Figure 3. The blue continuous line shows the mean density in 
equal volume regions as a function of the comoving distance when 
only the magnitude selection effect is taken in to account. The 
slight dependence of the density with distance shows evidence for 
a radial bias. The blue dot-slash-dot line is the linear fit adopted 
to model this radial bias and the red dashed line represents the 
mean density after correcting the density field by the magnitude 
and radial selection effects. 

cates that the effect here, is significantly less severe in our 
volume limited sample than in magnitude limited samples. 
To correct for this effect, we model the dependence of the 
mean density with comoving distance as S^dc) = a * d c + b 
with a = 0.0025 and b — 0.1565 for distances in Mpc. Tak- 
ing S2(d c ) into account, the radial trend in the mean value 
of the density field disappears, as shown in Figure [3] 

3.3 Boundary Effects 

Due to the large sizes of superclusters and to the limited 
volume of our galaxy sample, care should be taken to avoid 
boundary effects that can affect the analysis described in 
the next sections. 

With this aim, we considered only structures where all 
galaxies have comoving transversal distances from any vol- 
ume boundary border larger than a = 8 h _1 Mpc. Figure [3] 
shows the region occupied by our initial sample (in black) 
as well as its boundary points (in green). Superclusters in 
the region with 240° < a < 253° and -2° < S < +2° were 
excluded from our sample to avoid boundary effects. Ex- 
cluding superclusters at the boundary, our final superclus- 
ter sample has 880 structures above the threshold Dl and 
409 structures above D2. Their main properties are shown 
in Table 1. 



4 MORPHOLOGICAL DESCRIPTION OF 
SUPERCLUSTERS 

Since superclusters of galaxies are not virialized structures, 
they present a large variety of shapes. In this section we 
present the shapefinder method which we adopted to de- 
scribe the morphology of these structures. First we discuss 
the ellipsoidal model to describe the structures and then we 



Morphological Properties of Superclusters of Galaxies 5 




Figure 4. Region occupied by our sample on the sky (in black). 
Green dots represent the boundary of the region. 



use Minkowski Functionals to obtain a morphological pa- 
rameter. 



4.1 The ellipsoidal model 

A simple description of the tri-dimensional morphology of a 
body can be obtained through its best-fit tri-axial ellipsoid. 
This model has 8 free parameters: three for the centroid of 
the structure, two for its orientation and three for the semi- 
axes oi, CL2 and a-j,. 

The parameters of the ellipsoid can be inferred from the 
inertia tensor, i.e., the matrix of second-order moments of 
particle positions, 



^ L k W k XiXj, 

k 



(8) 



where x\ represents de i-th coordinate of the k th galaxy with 
respect to the object centroid. The matrix ly can be diag- 
onalized and the diagonal eleme nts are proportional to the 
best-fit ellipsoid semi-axes (e.g.. iJang-Condell fc HernquistJ 
l200ll ; IPlionis et af]|l99ll ; iKolokotronis et al.ll200 



h = 
h = 



(a 2 + a 3 ) 



— (ai +"3) 

EjLiW i f 2 2 
z (ai +a 2 ). 



(9) 



Solving the system of equations above, the three semi-axes 
are determined, with the assumption that a\ > 02 > a A . 

4.2 Minkowski Functionals 

Minkowski Functionals (MFs) represent an important tool 
to describe structures and objects, since they charac- 
terize their geometry. We follow here the formalism of 
ISahni. Sathvaprakash fc Shandarinl (119981 ) , which uses ellip- 
soidal models for the morphological description of the ob- 
jects. 

Having as input parameters the semi-axes ax, 02, 03 



obtained in the previous section, we can determine, for an 
object or isodensity contour, four parameters: the volume 
(V), the surface (S), the integrated mean curvature (C) and 
the integrated Gaussian curvature (Q), also called genus. 

The parametric equation for an ellipsoid with semi-axes 
01, a.2 and 03 can be written as 

r(8, <f>) = ai (sin#cosci)i + 02 (sin#sin<jf>) j + as(cos8)k. (10) 

We now define 
E — r e ■ re, 
F = r ■ t$, 
G = r<j,- r>, 
L = r e e ■ n, 
AI = r S • n, 
TV = r^ ■ n, 
where 

r-0 = dr/d<j>, 
re = dr/d6, 

= d 2 r/d(j> 2 , 
ree = d 2 r/dd 2 , 
re4> = d 2 r/d9d(f>. 

The vector n represents the unit vector perpendicular to the 
surface and is defined as 

n = r g x r^/j r g x j. 

The four geometrical quantities can then be written as 



\JEG - F 2 d6d(t>, 



C = 



Q= t- 



ki + ki 



AS, 



-1 
4tt 



kik2dS, 



V : 



-11010203. 



(11) 



(12) 



(13) 



(14) 



The principal curvatures of the ellipsoid are fei and fc 2 , 
and the product and sum of these quantities are 

EN + GL- 2FM 



ki + k2 = 



k\k2 = 



LN 



EG 
M 2 



F 2 



EG- F' 2 



(15) 



(16) 



Three parameters are introduced, Hi, H2 and H3, 
which have dimensions of length: H\=3V/ S, H2=S/C and 
Hs=C /4n . Combinations of these parameters provide two 
important shapefinders, K\ (planarity) and K2 (filamentar- 
ity), which can be expressed as 

Kl = **Z*± (17) 



and 



K 2 = 



H 2 + H! 



H3 — H2 
H3 + H2 



(18) 
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Figure 5. The histogram shows the distribution of K1/K2 for 
SDSS superclustcrs considering Dl (solid line) and D2 (dotted 
line). Both distributions present a median value around unity. 



The vector K = (Ki , K2 ) has an amplitude and direction 
which determine the shape of a certain 3D surface. An ideal 
pancake-like object presents one dimension which is much 
smaller than the others, so Hi « H2 — H3 and conse- 
quently K ~ (1,0). For an ideal filament, Hi >> Hi — Hz 
and so K ~ (0, 1). Considering surfaces like ribbons, the 
shapefinders have three distinct dimensions, i.e., Hi « 
H2 « Hz and K ~ (a, a) with a < 1. It is worth mention- 
ing that, for a sphere, Hi — H2 = Hz and hence K = (0, 0). 

This formalism can be used to classify objects with dif- 
ferent shapes, so we consider two morphologies: 

• objects with K1/K2 > 1 are classified as pancakes. 

• objects with < K1/K2 < 1 are classified as filaments. 

The range of K1/K2 for ribbons is somewhat arbitrary 
(K1/K2 — 1), so we decided to exclude this morphology 
from our classification. The shape statistics K1/K2, through 
the so-called "shape spectrum" , was first applied to astron- 
omy by B01. 

Figure [5] shows the distribution of the morphological 
parameter K\/K2- Table 1 presents some statistical prop- 
erties of the supercluster morphologies considering the two 
threshold densities discussed in Section 3.1. 



5 RESULTS AND DISCUSSION 

In this section, we investigate how overall properties of su- 
perclusters, like their richness and total luminosity, corre- 
late with their morphological properties. We also repeat the 
analysis for structures extracted from numerical simulations, 
with the main objective of comparing the supercluster prop- 
erties in the position and velocity spaces. 



5.1 SDSS superclusters 

Figure [S] and Table 1 indicate that we have essentially the 
same number of objects classified as filaments or pancakes 
in our sample, for the two density thresholds discussed here. 



Figure 6. Richness (top) and total luminosity (bottom) of 
SDSS superclusters as a function of the morphological parameter 
K1/K2 for the density threshold Dl. The median and quartiles 
of the distribution are shown. 



We found 436 filaments and 444 pancakes for the threshold 
Dl and 204 filaments and 212 pancakes for the threshold 
D2. Hence, our result, does not confirm previous works (e.g. 
B01), where a prevalence of objects classified as filaments 
was found. 

In order to investigate relations between morphology 
and properties of superclusters, the Spearman rank-order 
correlatio n coefficient, r s , w as used to measure possible cor- 
relations ijPress et al.ll2007h . We have also computed the 
two-sided significance level of the null hypothesis of ab- 
sence of correlation (or anti-correlation), P(Ho); a small 
value of P(Ho) is indicative of strong correlation (or anti- 
correlation) . 

Figure [6] shows the richness and total luminosity of 
SDSS superclusters as a function of the morphological pa- 
rameter K1/K2, considering the threshold Dl. The lines 
represent the median and quartiles for each bin of K1/K2. 
Despite the large scatter, there is a significant trend be- 
tween richness (and consequently total luminosity) and the 
morphological parameter Ki/Ki- In both cases we found 
r s ~ —0.25 and P(Ho) < 10 -4 , showing that filamentary 
structures tend to be richer and more luminous than pan- 
cakes. A similar behavior is found for the threshold D2, with 
a correlation coefficient r s ~ —0.25 and P(Ho) < 10~ 4 for 
both richness and total luminosity. These results indicate 
that the trends of richness and total luminosity with the 
morphological parameter are not strong but are statistically 
significant. 

We have also compared the luminosity distributions of 
filaments and pancakes through the Kolmogorov-Smirnov 
(K-S) test (Figure 0. We conclude that their distributions 
are statistically distinct, presenting a K-S probability lower 
than 10 -3 that the cumulative luminosity distributions of 
filaments and pancakes are drawn from the same distribu- 
tion. A similar result is achieved using the threshold D2, in 
this case with a probability < 10 -3 . Indeed, Figure [7] sug- 
gests that the luminosity distribution of filaments is signif- 
icantly broader than that of pancakes, resulting in a higher 
number of filaments at high luminosities, in agreement with 
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Table 1. Statistics for observed and simulated superclusters for the two threshold densities. The table shows the 
number of superclusters classified as filaments (Nf), pancakes (N p ), the mean number density (n,sc) of superclusters 
and median values for total luminosity, richness, and K1/K2. For simulated superclusters, we also present our results 
in position space. 
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Figure 7. Observed luminosity distribution of filaments (contin- 
uous red histogram) and pancakes (dotted blue histogram), for 
density threshold Dl. 



Figure [6] A similar trend for superclusters classified as fil- 
aments to be richer and consequently more luminous had 
already been reported by E07b. 

Considering now only the brightest structures, those 
with log(L/Z/0) > 12.5, we may notice the prevalence of 
filaments over pancakes: 60.4% and 39.6%, respectively, for 
the threshold Dl and 62.5% and 37.5% for the threshold D2. 

Figure [8] shows the distribution of the semi-major axis 
resulting from the ellipsoidal fitting to the structures. Fil- 
aments and pancakes have different sizes, as demonstrated 
also by Figure 1101 with filaments comprising most of the 
largest superclusters. 
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Figure 8. Distribution of the semi-major axis resulting from the 
ellipsoidal fitting to the structures. Left: filaments, right: pan- 
cakes. Dotted lines correspond to the distributions obtained from 
the superclusters extracted from the numerical simulations de- 
scribed in Section |5. 21 Results for the density threshold Dl. 



5.2 Analysis of simulated superclusters 

We have repeated the above analysis for the four mock SDSS 
catalogues described in Section 2, applying the same proce- 
dures described in Sections and U 

Table 1 shows the number of superclusters identified in 
the simulations, classified as filaments or pancakes, as well 
as the mean number density and median values of richness, 
total luminosity and K\ / K2 for each sample of superclusters 
and threshold densities. Comparing the median values of 
total luminosity and richness from both threshold densities, 
slightly higher values are found in D2 because only richer 
and more luminous superclusters are identified using this 
threshold, increasing the median values. 

The same behaviour of observed richness and total lumi- 
nosity with the morphological parameter K1/K2 is present 
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Figure 9. Richness (top) and total luminosity (bottom) of su- 
perclusters extracted from numerical simulations as a function 
of the morphological parameter K\ / K2 for the density threshold 
Dl. The median and quartiles of the distribution are shown. 
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Figure 10. Semi-major axis resulting from the ellipsoidal fitting 
to the structures as a function of the morphological parameter 
K1/K2 for the density threshold Dl. The median and quartiles 
of the distribution are shown. Top: SDSS sample. Bottom: su- 
perclusters extracted from the numerical simulations described 
in Section 15.21 



in the simulations. Figure [9] shows the relation between rich- 
ness and total luminosity and K1/K2 for the superclusters 
extracted from the simulations with density threshold Dl. 
The value of the Spearman coefficient is r 3 — —0.25 and 
P{H ) < 1(T 4 considering the threshold Dl, and r s = -0.25 
and P(Hq) < 10~ 4 for D2. This anti-correlation between 
richness and total luminosity and K1/K2 in the mock cat- 
alogues is very similar to that obtained for observed super- 
clusters (Figure [6]), indicating a good agreement between 
simulations and observations. 

Peculiar velocities may have an important influence 
on supercluster identification. The observed redshift of 
extragalactic objects has a component produced by the 
Hubble flow plus a peculiar velocity component due to 



f 9Olo , [ea >30% / 


f golo<le5 >60% / 


f 9oloile5 >90% / 


■ (r 5 ,P(H )) = (0.424,<10"^i 


■ (r s ,P(H ))=(0-532,<10"^T 


- ( r .> p ( H o)) = ( u - 6u1 >< 10 /i> 


• • Y 


" * * 




• "V / 




* * 1/ 






Us* ' 


'vuKV ; 










jfte.C. • 




w9t • 





01 23401 23401 234 
K,/K Jp „, K./K 2p „, K,/K 2p „ 



Figure 11. Comparison between the values of the morphological 
parameter of simulated superclusters identified in the velocity and 
position spaces. Each panel corresponds to different limits in the 
fraction of galaxies in common (/ 9 ) in structures identified in the 
velocity and position spaces. 



gravitational interactions. Consequently, there is a differ- 
ence between measured distances in velocity and posi- 
tion spaces which produces some features in the velocity 
space such as the fingers-of-god and the Kaise r effect (e.g., 
iBahcall. Soneira fc Burgettl 119861 : lKaiserlll987T L In order to 
study the influence of peculiar velocities on supercluster 
identification and mo rphology, we identifi ed superclusters in 
the mock samples of ICroton et al.l l)2006h using two sets of 
redshifts available in the simulations: the first one including 
peculiar velocities (velocity space) and the second without 
them (position space). The influence of peculiar velocities 
on the identification and morphological classification of su- 
perclusters can be studied by comparing the properties of 
simulated superclusters in both spaces. Table 1 also presents 
the main features of simulated supercluster samples identi- 
fied in velocity and position spaces. Qualitatively, results are 
similar for density thresholds Dl and D2. 

Due to peculiar velocities, galaxies associated to a struc- 
ture in position space may or may not be associated to the 
same structure identified in velocity space. In order to com- 
pare the morphologies of superclusters identified in both 
spaces, superclusters from position space were associated to 
velocity space ones according to the percentage of galax- 
ies (f g ) in common (30%, 60% and 90%). Figure ITT1 shows 
Ki/K 2veloc versus Ki/K 2pos for different percentages of 
galaxies in common. In all cases there is a significant correla- 
tion between the morphological parameter measured in both 
spaces, despite the large scatter, which increases as f g in- 
creases. Interestingly, there is no significant bias in the mor- 
phological parameter measured in velocity space compared 
with the values measured in position space. This is actually 
due to the high values of the smoothing parameter adopted 
here (see Section [3}. As another test of how morphology is 
distorted in redshift space, we have considered galaxies in 
superclusters identified in redshift space and compared the 
morphology they trace in both velocity and position spaces. 
Figure Q2] shows the distribution of Ki /K 2 in velocity and 
position-spaces in this case. The median values of K1/K2 
are similar (Ki/K2 veloc = 1.01 and ifi/iG^ = 1.04) in 
both spaces. A K-S test indicates that the two distributions 
are consistent with each other. Hence, this test also suggests 
that peculiar velocities are not important for our morphol- 
ogy measurements. 

We now consider the luminosity distribution of super- 
clusters identified in the velocity and position spaces in the 



Morphological Properties of Superclusters of Galaxies 9 




Figure 12. Distribution of the morphological parameter of super- 
clusters identified in velocity (black) and position (green) spaces. 
The median values of K1/K2 (top left) are similar and the K-S 
test shows that the distributions are consistent with each other. 



mock catalogues. Figure [131 shows this distribution for fila- 
ments and pancakes for the density threshold Dl. For each 
morphological class the distributions in both spaces are quite 
similar, as confirmed by the K-S test. It is worth mentioning 
that these distributions in velocity space resemble very much 
the observed distributions (Figure 9). This is confirmed by 
a K-S test comparison of the filament and pancake distri- 
butions with K-S probability equals to 0.85 and 0.42, re- 
spectively, for the null hypothesis that the data were drawn 
from the same distribution. Both probabilities show that the 
distributions are not statistically distinct. 

The trend of the semi-major axis (di) of the ellipsoidal 
fitting of observed and simulated superclusters with the mor- 
phological parameter shown in Figure [TO] provides an indica- 
tion that the largest structures tend to be filamentary. The 
trend is very similar for the observed and simulated struc- 
tures. 

Our results indicate that simulated and observed su- 
perclusters present similar properties. This is also the case 
for the luminosity distributions of observed and simulated 
superclusters, as can be seen in Figures [71 and [T3l However, 
E06 have found an absence of very luminous superclusters 
in simulations compared with observations. Indeed, a com- 
parison between the observed and mock samples shows that 
objects brighter than \og(L / Lq) > 12.5 are at least twice 
more frequent in the observations than in the mock samples. 

Additionally, by comparing the luminosity distribution 
of superclusters classified as filaments and pancakes, we have 
noticed that these classes have very different luminosity dis- 
tributions. Since filaments tend to be richer and more lu- 
minous than pancakes, it is fair to suggest that these two 
morphological classes represent different evolutive dynami- 
cal stages of the large scale structure, with pancakes possibly 
evolving to filaments. 



5.3 Statistical reliability of supercluster 
identification 

Our method of supercluster identification is strongly depen- 
dent on the threshold density and on the smoothing param- 
eter. Thus it is necessary to verify its robustness by checking 
what fraction of our detections is expected to be real or just 
statistical noise. 

To investigate the expected num ber of random clu mps, 
we follow the approach proposed by iBasilakos! l|2003h . We 
ran a large number (200) of Monte Carlo simulations where 
we randomized the equatorial coordinates of the galaxies, 
keeping their comoving distances and the sample bound- 
aries in order to preserve the selection function. Since this 
randomization destroys the SDSS clustering, the structures 
identified in each simulation will be due to statistical noise. 
The probability of identifying real superclusters in our sam- 
ple through this procedure is then 



P= 1 



Nsdss ' 



(19) 



where N ran d is the number of structures identified in the 
randomized samples and Nsdss is the number of structures 
identified in our original SDSS sample. A probability close 
to 1 means that the number of spurious objects produced 
by our supercluster identification method is small. 

This analysis was performed for the same threshold den- 
sities as in section [3TT1 Figure [T4l shows the resulting proba- 
bility distribution obtained with the threshold Dl, present- 
ing a median probability P around 85%. For the threshold 
D2, the median probability is higher than 99.8%. Figure [T5l 
compares the luminosity distribution of our Dl SDSS super- 
clusters with that obtained from the simulations. The result 
clearly indicates that most random structures have luminosi- 
ties significantly lower than those of the SDSS superclusters. 
This result is also noticed by comparing the median values 
of the SDSS and random luminosity distributions. 

The same morphological analysis as described in sec- 
tion 2] was applied to the structures identified in the Monte 
Carlo simulations. No trends were found between the to- 
tal luminosity or richness and the morphological parameter 
K1/K2, indicating that the results displayed in Figure [S] are 
indeed real and not produced by random fluctuations or by 
the selection function of the sample. 

Another test of robustness of the main trends identi- 
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Figure 14. Distribution of the probability of identifying real su- 
perclusters for the density threshold Dl. The median value (red 
dotted line) is around 85%. 
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Figure 15. Luminosity distributions of SDSS supcrclustcrs (con- 
tinuous line) and random superclusters (dotted line) for the den- 
sity threshold Dl. The median of logarithmic values (top left) 
show that random superclusters present significantly lower lumi- 
nosities than SDSS superclusters. 



fied in this work can be done by considering only structures 
with total luminosities larger than 10 12 Lq , since our Monte 
Carlo simulations indicate that the expected number of spu- 
rious structures in this case is only ~3.7%. All our previous 
results are confirmed. 



output of the Millenium Simulation (jCroton et al.|[2006l ). to 
confront theory and observations and to examine the role of 
some systematic effects. 

Superclusters were identified by the density field 
method with an Epanechnikov kernel, taking into account 
selection and boundary effects. By comparing the results ob- 
tained with this kernel with those obtained with a truncated 
Gaussian kernel, we show that the results are not strongly 
dependent on the kernel choice. The kernel smoothing pa- 
rameter and the dimension of the density field cell were cho- 
sen as cr=8 h _1 Mpc and i C ei=4 h _1 Mpc. Two threshold den- 
sities were chosen to evaluate their influence on superclus- 
ter identification and morphology. The first maximizes the 
number of structures (Dl) and the second is chosen by lim- 
iting the size of the largest superclusters to ~120 h _1 Mpc 
(D2). We found that, at least qualitatively, our results do 
not depend on the density threshold used. Each superclus- 
ter is characterized by its richness and total luminosity, as 
well as by a morphological parameter determined through 
Minkowski Functionals, which allows their classification as 
filaments or pancakes. 

Following E07b, we have found significant correlations 
of the morphological parameter K1/K2 with richness and 
total luminosity in both the observed and mock supercluster 
catalogues, indicating that filaments tend to be richer and 
consequently more luminous than pancakes. 

To evaluate the influence of peculiar velocities on super- 
cluster morphology, we have used mock catalogues to exam- 
ine structures identified in the velocity and position spaces. 
We conclude that peculiar velocities do not play a signif- 
icant role in our results, probably due to the large kernel 
smoothing length adopted in this work. 

We have found a trend between supercluster total lumi- 
nosity (or richness) and morphology, with filaments being in 
the mean more luminous than pancakes. A similar behaviour 
was found by analyzing mock catalogues. The analysis of 
Monte Carlo simulations of randomized galaxy distributions 
indicates that these trends are real and are not produced by 
random fluctuations or selection effects. 

Finally, we compared the luminosity and size distribu- 
tions of filaments and pancakes. We have found that they are 
significantly different, with filaments presenting a broader 
distribution of sizes and luminosities. Again, similar results 
were obtained by the analysis of the mock catalogues, show- 
ing that filaments and pancakes represent distinct morpho- 
logical classes of superclusters in the Universe. Also, since 
filaments tend to be richer, more luminous and larger than 
pancakes, it is plausible to think that pancakes evolve to- 
wards filaments. 
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6 SUMMARY AND CONCLUSIONS 

We have selected a volume-limited sample of galaxies from 
the SDSS in order to identify galaxy superclusters, deter- 
mine their morphologies and investigate some of their fea- 
tures. Our SDSS sample contains galaxies with M r < —21 
in the redshift range 0.04 < z < 0.155 and covers stripes 
10 to 37. We have also analyzed simulated light-cones based 
on a semi-analytic galaxy evolution model applied to the 
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APPENDIX A: THE INFLUENCE OF THE 
SMOOTHING KERNEL ON SUPERCLUSTER 
PROPERTIES 

The kernel in Equation 2 plays a fundamental role in the 
density field calculation and therefore it is important to ver- 
ify how sensitive our results are to the kernel choice. We 
address this point by comparing supercluster properties ob- 
tained with the Epanechnikov kernel with their properties 
obtained with a truncated Gaussian kernel, defined as 

K(r,a) = ±e^, (Al) 

with a cutoff at 3-cr. 

For this exercise the density field was calculated as 
described in Section 13.11 but with a smoothing parameter 
cr=2.4 h _1 Mpc for the truncated Gaussian model. This value 
leads to 1038 superclusters for the threshold Dl with lumi- 
nosities similar to those discussed in Section [5.11 as shown 
in Figure IAT1 A similar result is obtained by comparing the 
richness of superclusters obtained with the two kernels. 

In order to compare the structures obtained with these 
kernels, we have matched the two supercluster catalogues 
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Figure Al. Comparion between the total luminosity of su- 
perclusters identified using the Gaussian (log(L/ Lq)q) and 
Epanechnikov (log(L/ Lq)e) kernels with parameters described 
in the text. 
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Figure A2. Distribution of the morphological parameter K1/K2 
for superclusters identified using the Epanechnikov (dotted line) 
and Gaussian (continuous line) kernels. At the top right we show 
the median values of each distribution. 



by assuming that the distance between the object center- 
of-mass in both catalogues is lower than 8 h Mpc. In this 
way, we have identified 368 objects in common. 

In Figure [A"2l we compare the distributions of the mor- 
phological parameter K1/K2 obtained with the two kernels. 
The distributions have similar medians and the K-S test did 
not distinguish them. The same trends present in Figure O 
were found using the Gaussian kernel, with a correlation co- 
efficient r s ~ —0.18 and P(Hg) < 10 -3 for both luminosity 
and richness. Similar results are obtained with the threshold 
D2. We conclude that the trends described in the text are 
actually robust with respect to the kernel choice, as far as 
sensible kernels are used. 



